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We calculate the superconducting critical temperature Tc, the singlet pair function ^'+(1), and 
triplet pair function '^~(x) of superconductor/normal metal/ferromagnet (S/N/F) trilayers using 
the linearized Usadel equation near Tc.. The Green's function method developed by Fominov et al. 
for the S/F bilayers is extended to the S/N/F trilayer systems. The S of the trilayers is taken to be 
an s-wave singlet pairing superconductor, and the S/N and N/F interfaces are modeled in terms of 
the interface resistances parameterized, respectively, by and 'y^^ . We present the Tc, 'i'^{x), 
and ^'"(x) for typical 7jf^, J^'^ , and the exchange energy Ecx'. (a) For a small (large) 'j^^ , Tc of 
S/N/F trilayers, as djv is increased, increases (decreases) on the length scale of N coherence length 
^jv with a discontinuity at djv — due to a boundary condition mismatch, (b) Tc{dF) shows a non- 
monotonic behavior like S/F bilayers with a weakened shallow dip. (c) The odd frequency triplet 
component ^"(x), induced by Ecx and proximity effects, has a maximum near the N/F interface 
and decreases on the length scale i^cx in F. It also penetrates into N and S regions on the length scale 
and 5s, respectively. Based on these results we make comments on the experimental observation 
of the odd triplet components and the recent Tc measurements in Nb/Au/CoFe trilayer systems. 

PACS numbers: PACS: 74.45. +c, 74.78.Fk, 74.25.Ha, 74.20. Rp 



I. INTRODUCTION AND MOTIVATION 

When different states of matter are in contact, one 
affects the other over the coherence length ^, which is 
referred to as the proximity effects.— '^'^'^ For instance, 
when a superconductor (S) and a normal metal (N) are 
in contact as in an S/N bilayer, the singlet pair function 
^'^{x) ~ {(j)^{x)(j)i{x)), where 4><j{x) is the field operator 
of spin cr electron at the position x and ( ) means the 
thermodynamic average, penetrates into the N region and 
decays exponentially on the length scale of N coh erence 
length ^jv. For diffusive metals, = xAD/v/SttT where 
Djv is the diffusion constant of the normal metal and T 
is the temperatureJiiS Note that, however, the gap en- 
ergy given by A(a;) = V'^^{x) vanishes in N because the 
pairing interaction = in the N region. The supercon- 
ducting temperature Tc of an S/N bilayer as a function 
of N thickness cIm also decreases exponentially on the 
length scale of ^at. Here, Tc is defined with the current 
parallel to the interfaces. Therefore, Tc of an S/N bilayer 
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is determined by the highest temperature at which the 
singlet pair function '^'^{x) does not vanish at least at 
one point within the bilayer. 

The proximity effects in S/F bilayers have also been 
studied intensivelyj^'"* Like the S/N case, ^'+(2;) of an S/F 
bilayer decreases exponentially in the F region, and Tc de- 
creases as the F thickness dp increased. Many of the 
interesting S/F proximity effects occur at the nanoscale 
range of layer thickness. The observation of these ef- 
fects was made possible by the progress of fabrication 
technique of high quality hybrid S/F layers. The S/F 
proximity effects have the following differences compared 
with S/N: (a) The decay length scale in F, ^es, given by 
£,ex — y'Dp/ Eex , is much shorter than ^pf because the 
exchange energy E^x is much larger than T. (b) The pair 
function oscillates as well as decreases in the F region be- 
cause the net momentum of a Cooper pair is non-zero in 
Ff^ By the same reason Tc{dp) does not monotonically 
decrease like S/N as dp is increased but shows a shallow 
dip. The oscillation has been confirmed by the observa- 
tion of the "TT-state" in S/F/S Josephson junction,— '^'^'^ 
where two superconductors separated by a ferromagnet 
of appropriate thickness can have a phase difference of tt. 
(c) Triplet pairing component (TPC) is induced in addi- 
tion to the dominant singlet pairing component (SPC). 
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The triplet pairing is realized in ^He and Sr2Ru04, for in- 
stance, in the form of p- wave pairing. The induced TPC 
in S/F systems in the diffusive limit, however, is even (s- 
wave) in momentum but odd in frequency, referred to as 
odd frequency triplet pairing component (OFTPC). The 
OFTP was first suggested by Berezinskii in the context of 
•^He superfluidityjil but it turned out that it is the p-wave 
triplet pairing that is realized in ■^He. The OFTPC is in- 
duced from s-wave SPC because time reversal symmetry 
is broken by the exchange field in S /F bilayers. Just like 
non-s-wave pairing components are induced from s-wave 
one when the translational symmetry is broken, OFTPC 
is induced in S/F bilayers from the singlet component 
because the time reversal symmetry is broken by the ex- 
change field in F region. This can be seen directly from 
Fig. m where the induced TPC vanishes as Eex — > 0. 

While it is natural to expect that the characteristic 
length scale in F region is very short given by £,ex, it 
was realized that there exist long range proximity ef- 
fects in S/F systems set by the long length scale of 
= ^ Dp 12-kT . For instance, the long range triplet 
pairing was suggested to understand the unexpected con- 
ductance increase of S/F structures below Te J^^i^^'^'* The 
condition for the occurrence of the long length TPC is 
that the magnetization is not unidirectional. When the 
magnetization in F is unidirectional, the induced triplet 
pairing component has the short coherence length ^ea, like 
the singlet one and can be represented as S'z = com- 
ponent along the magnetization direction in Fj^^ More 
explicitly, it is given by (^t (x, t)0x (a;) (/>x(a:,r)0|(a;)), 
where r is the imaginary time. When the magnetization 
is not unidirectional, on the other hand, other triplet 
components with — ±1 given by ((/>| (x, t)(;/)| (x)) and 
(0x(x,t)(/)|(x)) are also induced because 5*2 = along a 
spin quantization axis is not pure Sz — ^ but a mixture of 
Sz = 0, -1-1, and —1 along a diff'erent axis. The = ±1 
TPC has the long coherence length (^p in F because the 
exchange field is not pair breaking for them. 

We here note the following distinctive features of the 
TPC which appears in S/F systemsi^ See below for more 
detailed discussion, (a) The S in an S/F bilayer is a 
conventional s-wave singlet superconductor. The exotic 
state of OFTPC is induced not because the S of the 
S/F is exotic. It is induced because of the proximity 
effects and the broken time reversal symmetry, (b) TPC 
is induced in addition to the dominant SPC, and may 
have short coherence length ^e^, or long coherence length 
^p. (c) TPC is even in momentum but odd in frequency, 
(d) The Tc of an S/F bilayer is determined by the SPC 
(See Eq. ^ below.), and TPC can change indirectly 
only by changing SPC through the boundary conditions. 
See Sec. II in the following for detailed discussions. The 
OFTP is an exotic state of matter, but it has not been 
directly observed yet. A lot of experimental and theoreti- 
cal work has focused on finding fingerprints of the OFTP 
correlationsiiSiii The basic idea is that one convert the 
short range triplet component into the long range one by 
controlling or inducing nonuniformity of the F magneti- 



zation direction, and observe the long length scale triplet 
component with appropriate experiments. Even when 
the exchange field is uniform in the F, the long range 
components are induced when the interface is spin ac- 
tive which causes spin flip or spin rotationJ^iiSii^iSi The 
present formulation can be straightforwardly applied in 
this case too with a modified version of the boundary 
conditions (BC) given by Eq. (|13p. Recently, there was a 
report which strongly hints the existence of long length 
OFTPC in S/F/S^ 

The superconductor-ferromagnet systems, as described 
above, exhibit very interesting and exotic behavior, and 
it will be important to understand them in a systematic 
way. To do that, we have theoretically studied the S/N/F 
trilayer systems. By varying the intervening N thickness 
djsi and the interface resistances between S and N and 
N and F, represented in terms of and 7^^, respec- 
tively, we can control the proximity effects in a system- 
atic way. In this paper, the proximity effects in S/N/F 
trilayers are studied for Tc, SPC, and TPC. This paper 
is organized as follows: After this section, we present in 
Sec. II the generalization to S/N/F trilayers of the nu- 
merically exact Green's function method developed by 
Fominov et alM for S/F bilayers. The resulting Usadel 
equations subject to appropriate boundary conditions are 
then solved self-consistently by numerical iterations for 
Tc, SPC, and TPC. The results wifl be presented in Sec. 
Ill for typical values of 7^^, 7^^''^, and E^x- Based on 
these results, we will make some comments on the exper- 
imental observation of the odd frequency triplet compo- 
nents and the recent Tc measurements in Nb/Au/CoFe 
trilayer systemstS^ We then conclude with the summary 
and concluding remarks in Sec. IV. In Appendix, we col- 
lect some technical details for calculating the singlet and 
triplet pairing components. 



II. USADEL FORMULATION 

A. Usadel equation for S/N/F trilayers 

We wish to understand the proximity effects in the 
S/N/F trilayer systems, which is schematically shown in 
Fig. HI The S is a conventional s-wave singlet pairing 
superconductor like Nb, N is Au, Cu, or Al (above Al's 
Tc), and F is CoFe or Ni. Each of S, N, F is characterized 
by the coherence length ^ and resistivity p. In addition, 
S layer is described by superconducting critical tempera- 
ture Tco, and F by the exchange energy E^x- The layers 
of thin films are in the dirty limit and we employ the 
dirty limit quasi-classical theory, the Usadel equation, to 
describe them;^ Moreover, since we are interested in cal- 
culating the superconducting transition temperature Tc 
and the pair functions near T w Tc of S/N/F, we will 
use the linearized Usadel equation. Near Tc, The normal 
Green function is G = sgn{ijjn)5a,a' , where a and a' are 
the spins of two electrons forming a Cooper pair, and the 
anomalous function is -F ^ 0. The Usadel equation is 
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therefore linearized with respect to F. In a general case 
where all four components of F are kept, F may be rep- 
resented as a 2 X 2 matrix or a 4 component vector. In 
the present problem where the exchange field is uniform 
over the F region, only two components, SPC and short 
range TPC, appear in the equation and the anomalous 
function F may be represented as a two component vec- 
tor or scalar complex function. Because we will consider 
long range TPC in a subsequent study where a vector 
representation is convenient we will formulate the present 
problem in terms of two component vector 

where i = S, N, or F, and x is the coordinate perpen- 
dicular to the interface and the translational symmetry 
is assumed parallel to the interface. The superscript ± 
refers to the even and odd functions of the frequency. 

(2) 

where the anomalous function F is defined by 



F{x,iuj„) 



dTe'^-^Tr{^^ix,T)4>^ix,0)), (3) 



where Tr is the t ordering operator. Then, F{x, —iLOn) = 
Jp'' dre'""'^rT-((/);(a;, r)^l (x, 0)), and the even function 
F+ represents the SPC and the odd function F^ the 
TPC of the anomalous function. We will also use the 
SPC (TPC) to stand for ^+{x) (^-(a;)) which is just 
the summation over Matsubara frequencies of F^{x, iun) 
{F~ {x,iujn)) as given in Eq. ([50|l . 




FIG. 1: Schematic drawing of an S/N/F trilayer. The coor- 
dinate X is taken as perpendicular to the interfaces. 

The linearized Usadel equation near T Tc takes the 
following form.i^ 

i'^s'^Tc^^'Fs{x,iu)n) = |a;„|Fs(x,iw„) - A(a;), 



^]yTTTc-^FNix,iuJn) = \uJn\F N {x , iuJn) , 



{-d-N < X <0), (5) 



S,pT:Tc-^FFix,iuJn) = |w„|F_f(x,zlj„) 



' dx^ 



-i sgn (w„) Ecx(TiFFix,iUn), 

{—dp — djv < X < — djv)- (6) 



where u;„ = 2TrT{n + ^) is the Matubara frequency and 



CTi is the Pauli spin matrix. A(x) = 



A(x) 




Eq. 



explicitly shows that the exchange field E^x couples the 
SPC F+ and TPC Fp and induces TPC from SPC. The 

coherence length is ^ = \J 2ttT dirty limit where 

D is the diffusion constant. Then, 



Ds 



i 



N 



D 



N 



27rT ' 



6 
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(7) 



in each S, N, and F layer, as discussed in Introduction. 
The singlet component Fg determines the gap function 
A(x) 



A{x) = Vg{eF)TiT, V F+{x,ioj^) 



E 

a)„>0 



(8) 



where V and gi^ep) are, respectively, the pairing interac- 
tion and the density of states per spin at the Fermi level. 
It satisfies the self-consistency relation 



E 



A{x) 



F^{x,iLL!n) 



(9) 



where Tc is the superconducting critical temperature of 
the S/N/F trilayer and Tco is that of the S layer alone 
without the N and F layers. Eq. ^ explicitly shows 
that Tc is determined by the singlet component alone. 
The triplet components may affect Tc only by changing 
the singlet component through the boundary conditions 
discussed below. 

Being a differential equation, the Usadel equation must 
be supplemented by boundary conditions. The appropri- 
ate BC for the Usadel equation are well established and 
given as follows 3^ 



dx 



Fp{-dF-dN)^0, 



dx 
dx 



^N) 



iNF^F-T-F f{- 

dx 



0, 



ISN^N-r^Ni^) 

dx 



FM{~dN)~Fp{-dN 
FAr(O) 

Fsids) = 0, 



0, 
d 



^-d-x^-^- 



Fs(0) 
d 



dx 



(0<x <ds), (4) where 



PNi 



N 



PF^F 



ISN 



PnS.n 



(10) 

(11) 
(12) 

djv),(13) 
(14) 
(15) 

(16) 



The 7^ ^ and are dimensionless quantities charac- 
terizing S/N and N/F interfaces, respectively, given by 



lb 



NF 



PF^F ' 



lb 



SN 



PN^N 



(17) 
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where R^^ 



and i?f ^ 



are the interface resistances be- where 
tween N and F and S and N, respectively, and A is the 
interface area. 7^ ^ and are modeled in terms of 

spin conserving potential barriers (2^*2^ and do not intro- 
duce any spin modifying mechanism. On the other hand, 
the boundary of the F can be spin dependent or spin 
active^>iai20^ and the BC of Eq. ^ should be modi- 
fied accordingly. The spin-active interface can induce the 
long range triplet components. One possible description 
of the spin-active interface is that the parameter is 
given by a 4 X 4 matrix instead of a scalar with the state 
of Eq. Il]) enlarged to a four component vector. All of the 
four components, two of which are the long range ones, 
are non- vanishing in general because they are coupled via where 
the BC Eq. 



Af = 



Af = 



RcAp 
-ilmAF 

kp^F tanh kpcLF 



—HthAf 
Re^F 



(22) 
(23) 



Then, utilizing the BC at a; = 
the FAr(-djv) and £Fjv(-d7> 
oi¥F{~dN) and ^YFi-d^) 
find 



-dN of Eq. ini) and (fT3|) . 
) can be written in terms 
Then, using Eq. (PT|) . we 



'dN), 



A 



NF 



1 



AFir) 



InfAf- 



(24) 



(25) 



B. Boundary condition in terms of Fs 

The Usadel equation (g]) , © , © together with the 
BC (Iini), in]), (ini), (113, (ini), (IISl) and self-conslstency 
equation ([9]) forms a complete set of equations to describe 
an S/N/F trilayer near T « T^- It is solved by extending 
the numerically exact Green's function technique devel- 
oped by Fominov et al. for the S/F bilayers32i22, To uti- 
lize the technique, everything should be written in terms 
of Fs alone in the S region. The Usadel equation and 
the self-consistency relation are already written in terms 
of Fs alone in the S region, and the remaining task is to 
write the BC at a; = and a; = ds in terms of Fs only. 
The BC at a; = ds is already given by Eq. ([TS]), and 
the one at x = can be derived by writing the Usadel 
equation with the corresponding BC successively starting 
from the F to N regions. The procedure is very similar 
to Fominov et al^ 

First, for the F region we solve the homogeneous equa- 
tion dHD with BC (dni) to obtain 



_ cosh kF{xNF) cosh kp{xNF) 
^ ' -coshkpixNF) coshk*p{xis[F) 



Second, the anomalous function F^r in the N region 
described by Eq. ([5]) may be written as 



where 



xnf = X + d 



N 



nTr 



(19) 



Fjv(x, iujn) ~ F7v(— d^v) coshfcAr(a; -I- djy) 
, ^ 17 / J sinhfcAr(a; + djv) 

+^N-^t N(~dN) ; T , 

ax k]srt,N 



where 



N 



Then, we use Eqs. and to obtain 



^iV:^Fjv(0) = ^JVsFAr(0), 

ax 



where 



Ans = ( 



= {An + Anf 1 



AnfAn 



(26) 



(27) 



(28) 



(29) 



An = kj\j£_j\j ta,iihkjs[dN- (30) 

Third, using BC at a; = given by Eqs. ([12]) and (fT4| . 
and Eq. (|28p . we can finally obtain the BC in terms of 
Fs(0) and £Fs(0) such as 



dFs(O) 
dx 



^sFs(O), 



(31) 



Note that the kp implies another length scale S,ex much 
shorter than ^f alluded earlier. Because E^x ^ |w„| for 
relevant Matsubara frequencies, one has 



where 



-4s = 1 



An sir)' 



IsnA 



NS- 



(32) 



6 



1 



Re{kF) 



Df_ 



(20) Eq. ([3T]) connects 



The length scale of pair function is given by ^ex in F 
region. From p8)) we obtain 



dF_ + 



(0) 



dx 



and 



rfFg (0) 
dx 



with 



F+(0) and 



Eg (0). But, we need the BC which connect 



dF+{0) 



with 



Eg(0) only to solve the Usadel equation in S region for 
< a; < ds given by Eq. (g]). 



^F-j^FFi-dN) = ApFFi-dN), 



(21) 



52 



\Un\F+{x,lLOn)^~A{x). (33) 
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We write from BC of 

Fg{x, iujn) = Cg{iujn) coshksix - dg), (34) 

where 



The solution of (l33l) is then 



ds 



F+{x,iujn)^ dyGix,y;iLU„)A{y). (47) 



ks 



V 



(35) 



Using this to ehminate the Fg (0) from ([HT|) for — 
we obtain 



dx 



(36) 



The gap function A(x) is determined by Eq. ([8]). 

Substituting (|i7|l into ([5]) gives the self-consistency 
equation, 

A{x)=Vg{eF)7TnJ2 [ ' dyG{x,y;iu;^)A{y). (48) 

From this self-consistency equation, we can write down 
the equation for Tc with respect to Tcq. 



where 

W{iuJn) = 757V 



A(a;) In 



Tco 



Asjj^'^ + ReBsN)+7SN 



ds 



dy M{x,y)A{y), 



(37) 



and As and Bsn are defined by 



tanh kf^xo 



As = ^5^5 tanh fcsds, (38) 
Bsn = [fcAT^AT tanhfcjv (c?Af + xo)]""^ (39) 

1 iNF 



kN^N + A-/ ■ 



(40) 



Green's function method 



Now, the problem is reduced to solving Eq. ([55)1 with 
BC ^ and To do it by the Green's func- 

tion technique, one needs to solve the following source 
equation.— 

nT,es-^Gix, y) - |c.„|G(x, y) = ^5{x - y), (41) 
with the BC corresponding to and 

^s^G{0,y) = W{iu:n)G{0,y), (42) 



is^G{ds,y) = 0. 



(43) 



The Green's function can be constructed as foUowsi^ 

ks/\u!n\ 



G{x,y]iuj„) 



sinh ksds + ^^^'^f cosh kgds 

x< ^i;^)^2(y), 0<x<y 

vi(y)V2{x), y < X < ds ' ' 



where 



vi{x) — cosh ksx-\ ^ "'^ sinh ksx, (45) 



ks^S 

V2{x) = cosh ks {x — ds)- 



(46) 



LO„>0 



6{x - y) 



G{x,y;iujn) 



(49) 



The above integral equation can be reduced to an eigen- 
value problem after we change the integration over y into 
a summation by discretization. The Tc is determined by 
the smallest eigenvalue of the discretized M{x,y) matrix, 
which yields the highest Tc. One plugs an arbitrary initial 
Tc into the right hand side of Eq. (|49p and calculates a 
new Tc, and iterates until the input and output are equal 
within a tolerance. After a convergence is reached, Tc is 
obtained, and the corresponding eigenvector to the small- 
est eigenvalue gives A(a;). Then, the anomalous function 
F+ can be obtained from Eq. (I47|) . and the singlet and 
triplet pair functions "^^{x) and \E'~(x) can be calculated 
from 



(50) 



Technical details for calculating SPC and TPC are col- 
lected in Appendix. 



III. RESULTS: AND PAIR FUNCTIONS OF 
S/N/F TRILAYERS 

A. Perfect interfaces 

Tc of S/N/F trilayers were calculated by self- 
consistently solving Eq. by numerical iterations. 

And Tc of S/N bilayers were calculated for comparison. 
The singlet and triplet pair functions were also calcu- 
lated as detailed in Appendix. Let us first consider the 



ideal case of perfect interfaces with 7^ = 7b =0- We 
take ps = 15.9627 /ificm, pN = 2.0 /zficm, pp = 40.0 
/ificm, ^5 — 7.0 nm, ^at = 110 nm, ^i? = 10.241 nm, 
Ecx = 1235 K, and Tco = 7.927 K, which are appropriate 
for S = 23 nm Nb, N = Au, and F = CoFe as reported 
in the reference22iS£. Recall that the condition EcxT ^ 1 
should be satisfied for the diffusive Usadel equation to be 
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FIG. 2: (a) Numerical results of Tc/Tco of S/N/F and S/N as 
a function of cLn, shown by dashed and solid curves, respec- 
tively. Both curves have the coherence length ^jv and merge 
as djv — > cx) as expected, (b) Tc/Tco of S/N/F and S/F as a 
function of dp, shown by the dashed and solid curves. Note 
the different length scales ^jv and ^ex between (a) and (b). 

applicable. Using D = ^VpT, and the i^^^ and Tco values, 
we obtain E^xT ~ 10^^ so that the Usadel equation is 
applicable to Nb/CoFe systems. 

We show in Fig. [UJa) Tc as a function of N layer thick- 
ness, dN, for ds — 23 nm and dp ~ 23 nm S/N and 
S/N/F. For S/N bilayers Tc(djv) decreases as djy is in- 
creased, but for S/N/F trilayers Tc{dN) increases as dN is 
increased. The length scale of both decrease and increase 
is given by the N coherence length ^n. As ^ oo the 
two Tc approach each other as expected. In Fig. ^h) we 
show Tc as a function of dp for S/F and S/N/F. Note 
the different length scales between (a) and (b) given by 
and S^ex, respectively. Tc{dp) for S/F shows the well 
known behavior of an exponential decrease with a shal- 
low dip. For S/N/F, we took dM — 50 nm. The presence 
of the N layer weakens the effects of exchange energy of 
F. Consequently, Tc{dp) for S/N/F, compared with S/F, 
is lower at c?i? = but higher as dp — > oo, and has a 
shallower dip. 

We may also calculate the pair functions as well as 
Tc from the Usadel equation for various S-N-F layers. In 
Fig. Ha) we show the SPC ^+{x) for S/N and S/N/F for 
ds — 23 nm, dN = 50 nm, and dp = 23 nm, normalized 



FIG. 3; (a) SPC of S/N/F and S/N shown by the solid and 
dashed lines, respectively, ^'(a;) is normalized by the value 
at the outer surface of the S layer, (b) SPC and TPC of an 
S/N/F trilayer shown by solid and dashed lines, respectively. 
TPC is magnified by 10 times for clarity. 

by the value at a; = 23 nm. It decreases monotonically 
in N but non-monotonically in F. The non-monotonicity 
of ^+(0;) and Tc{dp) results from the same physics. It 
may be understood as originating from the non-zero net 
momentum of a Cooper pair.'^ In (b) we show both SPC 
and TPC for an S/N/F. While SPC is maximum at the 
S outer surface of a; = 23 nm, TPC is maximum near the 
N/F interface. The position of TPC peak is roughly S^^x 
from the N/F interface. 

To see the TPC more clearly, we plot in Fig.|3]the SPC 
and TCP of S/F bilayers for several values of Eex- Note 
that for Ecx = F becomes like an N; the length scale £,ex 
becomes £,p and TPC vanishes. The TPC becomes small 
also for large Egx because TPC decays rapidly for small 
^ex oc Eex^"^- The TPC ^ as E^x ^ or ^ cxd, and is 
fully induced for an intermediate Ecx , as can be seen from 
figure (b). For each Ecx, the peak of TPC is determined 
by the competition between the exponential decrease of 
length scale ^cx and linear increase determined by the 
BC 

is^^si^) - lSF^F^Fp{0) = 0. (51) 
ax dx 

The TPC, therefore, has a peak at around ^cx away from 
S/F interface, which can be seen from Fig. 21 
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FIG. 4: SPG and TPC of S/F bilayers for several Ee^. The 
peak position of TPG is given by the decay length scale ^ex of 
SPG. Note that the TPG vanishes as F becomes an N {Eex = 
0. 



FIG. 5: Tc(dN) of S/N/F and S/N for several different inter- 
face resistances, (a) is for several with 7^^^ — 0. Two 
rc(d]v) curves of S/N/F and S/N merge as djv — » oo. (b) is 
for several 7^ ^ with 7^ ^ = 0. Note that for small 7^^^^ Tc 
increases as d]v is increased but decreases for large 76^^. 



B. Effects of interface resistances 

We now turn to the more realistic case with non- 
zero interface resistances. There are two interfaces in 
an S/N/F, between S and N, and N and F, parameter- 
ized by 7^''^ and 7^^, respectively. The main effect of 
the interface resistances is to weaken the proximity ef- 
fects and increase the Tc of S/N/F trilayers. This can be 
seen from Fig. [5j Figure (a) shows Tc vs. (In of S/N/F 
for several 7^^ with 7^^ = 0. Tc vs. cIn of S/N with 
the same 7^ ^ are also shown for comparison. Figure (b) 
shows Tc vs. djv for several 7^^ with 7^^ = 0. In both 
figures (a) and (b), Tc of a given increases as 7;, is 
increased. 

Another, perhaps more technical, way of saying this 
is from BC of Eqs. ([m, ([l2l), ([131), and (Ull). Non-zero 
7^''^ (7^^) introduces the discontinuity of the anomalous 
function at the interface. The anomalous function 
F+ and pair function therefore, increases in S but 
decreases in F region as shown in Figs. [6] and [T] An 
increase of \E'"'" (x) in S means an increase of Tc because it 
is the largest value of ^"'"(x) that determines Tc- (In the 
figures, however, ^'"'"(a;) is normalized at a; = ^5 and this 



point is not manifested clearly.) We plot in Fig.[n]^'^(a;) 
for S/N/F for several 7^^ values with 7^''^ = 0. And 
in Fig. [7] ^'^(a;) for several 7^^ values with 7^^ = 0. 
Note that '^^(x) in S region is increased as 7^ or 7^''^ 
is increased. 

Now, let us look at in more detail how Tc behaves as a 
function of for non-zero 7^^ or 7^^. When 7^^ = 
and 7^^ is increased for S/N/F, as shown in Fig. [SJ^a), 
Tc always increases as cIn is increased. However, when 
7^^ is increased with 7^^ = 0, Tc{dN) as a function 
of djv behaves somewhat differently. For instance, let 
us compare the two Tc{dM) curves for 7^^ = and 
Ib^ = 1.0. Tc{dN) increases with d^ for 7^-^ = but 
decreases for 7^^ = 1.0, and the two curves merge to- 
gether as diq — > 00. This is simple to understand. When 
7^^ = 0.0, the pair breaking effects of the exchange field 
of F fully influences S for small d^ , and its effect is weak- 
ened as c?Ar is increased. Tc, therefore, increases as d^ 
is increased. When 7^''^ = 1.0, on the other hand, the 
effects of F is almost blocked out and it behaves like S/N. 
Consequently, Tc of 7^^ = 1.0 is much larger compared 
with that of 7^''^ = 0.0 as d^ and decreases as d^ is 



8 



increased. When d^v ^ 200 nm, the effects of F are also 
blocked out irrespective of , and all curves merge 
together. 




-60 -40 -20 20 
X (nm) 



6: The SPC and TPC of S/N/F for several with 
= 0. The inset shows TPC in detail in F region. Non- 
produces a discontinuity in SPC and TPC at the 
S/N interface. 



FIG. 

zero 7f ^ 



We now turn to the calculations of pair functions for 
S/N/F. In Fig. H we took 7^^ = and plot SPC and 
TPC for several 7^ ^. Non-zero 7^^ introduces the dis- 
continuity in both SPC and TPC. The inset shows the 
peak and dip structure of TPC more closely. As dis- 
cussed for Fig. [31 the TPC has a peak around S^ex away 
from the N/F interface. In Fig.[7]we plot SPC and TPC 
for several 7^^^ with 7^^ — 0. The discontinuity is at 
the N/F interface in this case. 



Tc discontinuity due to boundary condition 
mismatch 



An interesting point for Tc{dM) of S/N/F trilayers ap- 
pears when we consider limit. Tc of this limit, of 
course, ought to coincide with that of the corresponding 
S/F of djv = 0. For example, of Nb/Au/CoFe tri- 
layers must approach that of Nb/CoFe as Au thickness 
goes to 0. Also, the change of Tc with djv should have 




-40 -20 

X (nm) 



FIG. 7: The SPC and TPC of S/N/F for several 7^ ^ with 
7^^ = 0. The inset shows the TPC in detail in F region. Non- 
zero 7^^^ introduces the discontinuity at the N/F interface. 



The TPC in F region is suppressed by 7;, 
F region or TPC in N region. 



less than SPC in 



the length scale oi £,n simply because there is no other 
length scale in the N region. We therefore expect that 
Tc{dN) of Nb/Au/CoFe trilayers should behave hke the 
curve shown by the long dashed line beginning with the 
Tc of Nb/CoFe represented by the solid square at the low 
left corner of the main plot in Fig. [51 The experimental 
Tc{df^) measurements, however, show quite different be- 
havior as shown by the solid squares reproduced in Fig. 
[8]^?. There is an abrupt increase of Tc with the length 
scale of approximately 2 nm as c? at is increased as shown 
in the inset, and Tc increases with the expected length 
scale as d^ is further increased. 

One may understand this experimental Tc vs. dif be- 
havior as follows. Recall that in the present formalism 
the interfaces are modeled in terms of the two parame- 
ters 7^^ and 7^^ irrespective of N thickness. S/N/F 
trilayers, as dN 0, still have two interfaces character- 
ized by 7^^ and 7^^^, but the corresponding S/F bilayer 
has one interface characterized by 7^ Unless the three 
7hS satisfy a special match condition, Tc{dN = 0) of the 
bilayer and lim^jv^o Tclf^iv) of the trilayers need not be 
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the same. This condition may be derived from the BC 
of Eqs. (HI]), Uni), and 1^. We add ^ and and 
take dN ^ 0. We then obtain using (fTTjl 



Tc{dN) as dj^ is increased was also observed in the epi- 
taxial Nb(110)/Au(lll)/Fe(110) trilayers by Yamazaki 
and his coworkers 4^ 



Fs(0)-F^^(0) 



7nf)^f^Fp{0). 



(52) 



For this to agree with the corresponding S/F bilayer, it 
is required that 



7^ 



NF 

■lb 



-I^INF. 



(53) 



This condition need not be satisfied. _l hc 

NF SN T c 

" or 7^ . In case of 



The may take 

an arbitrary value irrespective of 7^ 
a mismatch, Tc{dN) then shows a discontinuity at d^r = 
0. The right hand side is usually larger as is the case with 
the Nb/Au/CoFe trilayers, and of dN ^ S/N/F 
should be larger than that of S/F. 



: Rij^ Nb{23nm)/Au 


: Nb{23nm)/Au/CoFe(10nm)' 
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d„ (nm) 



200 



FIG. 8: rc(djv) of Nb/Au/CoFe trilayers and Nb/Au bilay- 
ers, represented by the solid and empty squares, respectively. 
Nb and CoFe layers are fixed at ds = 23 and dp — 10 nm. 
The inset is the detailed plot of the short length scale. The 
long dashed curve in the main plot is the expected behavior 
of Tc{dN) without the BC mismatch. The dotted line con- 
necting empty squares is a guide to eyes. Reproduced from 
the reference |23. 

Note, however, that in experiments d^ — > has a sin- 
gle interface of course, and Tc of Nb/Au/CoFe as dN ^ 
must agree with that of Nb/CoFe. As dN is increased 
from in Nb/Au/CoFe, there begin to form two distinct 
interfaces. For the electrons to feel the two interfaces as 
the theoretical model assumed, a finite width of N re- 
gion is necessary. This width of about 2 nm is the length 
which is needed to interpolate the Tc of Nb/CoFe and 
that of djv — * Nb/Au/CoFe. The theoretical step func- 
tion behavior of Tc{dN) then should appear as a short 
length scale of « 2 nm observed in the experimentsi^^ 
This length scale is not the material width of the Nb/Au 
interface which is approximately an order of magnitude 
smaller. The same short length of « 2 nm has also been 
observed for N = Cui^ Al^ The sudden increase of 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In this paper, we have extended the Green's function 
method developed for S/F bilayers by Fominov et al. to 
S/N/F trilayer systems. We have calculated the super- 
conducting transition temperature and the singlet and 
triplet pairing components near T k, Tc for S/N/F tri- 
layers. The interface between S and N (N and F) was 
modeled in terms of 7^^ {ib^)- They represent spin in- 
dependent interface resistances. As have been and will be 
reported separately, the experimental Tc measurements 
could be fit by the present formalism i^^i^^i^^i'^^ 

There seems, however, a room to improve the present 
formalism. One possible route might be to model the 
N/F interface with two spin dependent parameters as 
was done for example in ref. [13, or to model it in terms 
of a spin-active interfacej ^^'^^i^°'^^ The latter route will 
be an interesting problem in connection with the long 
range odd frequency triplet pairing components because 
they are induced even when the exchange field in F region 
is uniform if the N/F interface in spin-active. 

Also noteworthy in the S/N/F system is the observa- 
tion of the length scale of 20 nm in the Tc vs. djv 
measurements in Nb/Au/CoFe trilayers. The natural 
length scale in N region is the coherence length « 100 
nm as shown, for instance, in Figs. [5] and [5] In addition, 
the short length of ~ 2 nm was observed because of the 
boundary condition mismatch. In the same experiments 
the Tc as a function of dN also exhibited a small oscil- 
lation with the length scale of k, 20 nm. The origin of 
this is not clear. It seems difficult to understand this 
oscillation within the current Usadel formalism because 
there is only one length scale of S^n in the theory. It 
is interesting that this Tc oscillation was also observed 
in the epitaxial Nb(110)/Au(lll)/Fe(110) trilayers.^^ In 
the epit axially grown Nb(110)/Au(lll)/Fe(110) trilayers 
the Tc oscillation length scale was w 2 nm, which is an 
order of magnitude smaller than that of the sputtering 
grown Nb/Au/CoFe trilayers. More experiments with 
different conditions will be necessary to clarify the origin 
of the Tc oscillations as a function oi dN- 

We have calculated the singlet pairing component and 
short range odd frequency triplet pairing components in 
the present paper. The triplet pairing component is in- 
duced on top of the dominant singlet component by the 
proximity effects and time reversal symmetry breaking 
exchange field. In the present work, the exchange field 
is unidirectional over the F region, and only the short 
range triplet component was induced. If the orientation 
of the exchange field changes in the F region, or if the 
inner-interface of F is spin active, then the other long 
range triplet components will also appear. In the latter 



10 



case the procedures developed in the present manuscript 
can be directly applied to understand the long range odd 
frequency triplet components. We are currently inves- 
tigating the spectroscopic properties of them in various 
trilayers of superconductor and ferromagnet where the 
present formalism can be applied. 
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Note that A(a;) is the eigenvector corresponding to the 
smallest eigenvector M{x,y) of Eq. and Fg{x,iujn) 
can be obtained from Eq. (|47|) . We can therefore find 
Cg (iLUn) from Eq. (j59p . and calculate Fg {s,iLUn) using 
Eq. (|34p . This completes the calculation of the SPG and 
TPC in the S region. 

We now turn to calculation of the pairing functions in 
the N region. The function can be written as 

F^{x,iijJn) = F^(0,ia;„)coshfcjva; 



HN-^F^{Q,iuj,,) 



sinh fcjva; 
£.Nk_ 



N 



From the BC (dH, we get 



F^{x)=Fi{0) 



(1 



SN 



1 



As) cosh k^x 
sinh kj^x 



As 



£,NkN 



(60) 



(61) 



V. APPENDIX: CALCULATION OF SINGLET 
AND TRIPLET PAIR FUNCTIONS 

We will in this Appendix present the procedure to cal- 
culate the singlet and triplet pair functions using the Us- 
adel formahsm developed in Sec. H. As discussed there, 
the anomalous function Fs{x, iujn) in the S region can be 
written as 



F^{x,iujn) 
Fg{x,iuJn) 



F+{x,iujn) 
Cg {iujn ) cosh ks{x — ds) 



(54) 



The BC at a; = in terms of Fg alone is given by Eq. 



We also have 



where 



ax 



W{iuj^) 
-As 



(55) 



(56) 



(57) 



where W{u)n) and As are given by Eq. (f37|) . We therefore 
write 



{As - A) 



F+{0,iujn) 
Cg (iun) cosh ksds 



0. 



(58) 



Dividing this by Fg , we obtain 



" 1 ■ 







-i 



-As 



A. 





coshksds 



1 

-ks£,s sinh ksds 



W{iu!n) 
Cg{iuJn)/Fg{0, iUn) 



(59) 



Therefore we get 



N) 



-A, 



(1 



7^ 



ISN 

sinh kpfdN 



As) cosh k^d 



N 



ISN 



N 



FiiO) 

CiV^iv(l 



ISN 



-As cosh kpfd 



N 



SN 



ISN 



-A5) sinhfe^vrfw ^62) 



-F'^(— c?Ar) and ■^F^{—dN) is expressed as ^"5(0) and 
the results of equation (|59p is used together, we can get 
F±(x). 

Next, we calculate Fp{x) in F. The function Fp{x) is 
written as Eq. (HH]). We then use the BC (HI]) to write 



7nf£,f {kpCp sinhkpdp + kpCp sinh kpdp^ (63) 
d_ 
dx 

7nf£,f [—kpCp sinhkpdp + kpCp sinh kpdptj&A) 



We subtract and add (l63t and (l64l) to obtain 



+ _ ^N-^F^{-dN) " £.Nj^Fj^ i~dN) 



r+ — 



2^NP^pkp sinhkpdp 

^^^■^Jv(^'^Af) + ^N-^F^{-dN) 



2jNp£_pkp sinh kpdp 
We then obtain Fp by plugging Cp into Eq. (flSj) . 

^N-^F^j-dN) - S.N-£F^{-dN) 



(65) 
(66) 



Fp{x, iujn) = Re 



2"/NF^Fkp tanh kpdp 
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(cosh kpix + (In) + smhkpix + (In) tanh ^f^f) , (67) trilayers. The pair functions *(a;) are given by 



Fp {x,iajn) 



-Im 



^N£F+{-dN)-^N£F^{-d. 



^'jNF^FkF tanh fcirC^F 



^^{x)='KTe F^{x,iLOn). 



(69) 



(cosh A;F(a; + (In) + sinhfcF(a; + cIn) ta.nh.kFdF) ■ (68) This completes the calculation of the singlet and triplet 

pairing components over the entire region of an S/N/F 
We can therefore calculate the anomalous functions F in trilayer. 
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